Global classification of two-component 
approximately integrable evolution equations 



Peter H. van der Kamp 



Department of Mathematics, 
La Trobe University, Victoria, 3086, Australia 

Email: P.vanderkamp@latrobe.edu.au 



Abstract 



We globally classify two-component evolution equations, with homo- 
geneous diagonal linear part, admitting infinitely many approximate sym- 
metries. Important ingredients are the symbolic calculus of Gel'fand and 
Dikii, the Skolem-Mahler-Lech theorem, results on diophantine equations 
in roots of unity by F. Beukers, and an algorithm of C.J. Smyth. 

1 Introduction 

A long standing open problem is the classification, up to linear transformations, 
of two-component integrable equations 



where F^G are purely nonlinear polynomials in variables Ui,Vi, which denote 
the i-th x-derivatives of u{x,t),v(x,t). Among the many different approaches 
to recognition and classification of integrable equations, the so called symme- 
try approach has proven to be particularly successful, see for example [25j HQ] 
and references in there. Until recently, all results obtained were for classes of 
equations at fixed (low) order n. This situation changed dramatically when, by 
using a symbolic calculus and results from number theory, Sanders and Wang 
classified scalar evolution equations with respect to symmetries globally, that 
is, where the order n can be arbitrarily high 34J. Our aim is to obtain a similar 
result for the class of multi-component equations ([T]) . 

In the symmetry approach the existence of infinitely many generalized sym- 
metries is taken as the definition of integrability. A generalized symmetry of 
equation ([TJ is a pair of differential polynomials S = {Si, S2) such that equation 
^ is also satisfied by u = u + eSi, v — v + £82 up to order e^. This leads to 
the notion of Lie-derivative: C{K)S = <^ 5 is a symmetry of {ut,Vt) = K. 




(1) 



The Lie algebra of pairs of differential polynomials is a graded algebra. The 
linear part (aM„, bvn) has total grading 0, the quadratic terms have total grading 
1, and so on. Gradings are used to divide the condition for the existence of a 
symmetry into a number of simpler conditions: C{K)S = modulo quadratic 
terms, C{K)S = modulo cubic terms, and so on. This has been called the 
perturbative symmetry approach [5T]. In the same spirit the idea of an ap- 
proximate symmetry was defined |23j . If C{K)S = modulo cubic terms, we 
say that S is an approximate symmetry of degree 2. And, we call an equation 
approximately integrable if it has infinitely many approximate symmetries. 

We contribute to the above mentioned problem by globally classifying equa- 
tions ([T|) that are approximately integrable of degree 2. This is achieved by 
applying the techniques developed in the special case of so called S-equations, 
where any approximate symmetry of degree 2 is a genuine symmetry |12j . It 
extends older results obtained by Beukers, Sanders and Wang [21 [3]. See [17] for 
an overview on the application of number theory in the analysis of integrable 
evolution equations and |24j for more recent results. The present article is a 
revised and extended version of the report |13j . 

As remarked in |23j the requirement of the existence of approximate sym- 
metries of degree 2 is very restrictive and highly non-trivial. On the other hand 
an equation may have infinitely many approximate symmetries of degree 2, but 
fail to have any symmetries. This problem involves conditions of higher grading 
and is left open. 

2 Generalized symmetries 

A symmetry-group transforms one solution to an equation to another solution 
of the same equation. We refer to the book of Olver [25] for an introduction to 
the subject, numerous examples and applications. 

We denote A = C[u, w, -ui, wi, . . .] and q = A®A. We will endow g with the 
structure of a Lie algebra. For any K = {Ki^K^) & g the pair S = {Si, S2) € g 
is a generalized symmetry of the two-component evolution equation 



vanishes. Here 6q is the prolongation of the evolutionary vector field with 
characteristic Q, cf. [2^ equation 5.6], 




(2) 



if the Lie derivative of S with respect to K, 
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and the total derivative 

^ d d 

k—O 

The Lie derivative is a representation of g. This property, with P,Q(zQ, 

C{C{P)Q) = C{P)C{Q) - C{Q)C{P) (4) 

corresponds to the Jacobi identity for the Lie bracket [P, Q] — C{P)Q which is 
clearly bilinear and antisymmetric, cf. [29, Proposition 5.15]. Another way of 
expressing Q is saying that g is a g- module. Another g-module is given by A, 
the representation being C{K)F — 5k{F) with K ^ ^ A. 

The word 'generalized' stresses the fact that the order of a symmetry can 
be bigger than one. Generally symmetries come in hierarchies with periodic 
gaps between their orders. For example, the Korteweg-De Vries equation ut = 
U3 + itui possesses odd order symmetries only. Concurrently, the KDV equation 
has approximately symmetries at any order. 



3 Grading 

Denote cr„ = (u, 0) and = (0, v). HP in some g-module is an eigenvector of 
'C(cru) (or of £(tTt,)), the corresponding eigenvalue is called the u- (or v-) grading 
of P. If P has u-grading i and u-grading j we say that « + j is the total grading 
of P. One verifies that g can be written as the direct sum 

t>0 -l<i<t+l 

where elements of g*'-' have w-grading i and u-grading j. For example, {uiV2, V3V4) € 
g"'^ has total grading 1. Similarly we have 

A — 4* A* — Ai,i-i 

t>0 Q<i<t 

The crucial property of a graded Lie algebra is that the u-, (or v-, or total) 
grading of C{P)Q is the sum of the u-, (or v-, or total) gradings of P and 
Q. This follows directly from equation Q. Gradings are used to divide the 
condition for the existence of a symmetry into a number of simpler conditions. 

We study evolution equations of the form 

A'" + + • • • 

A'"^" + K-^^^ + A'"^i + K^^° + if 2.-1 + ... ^ (5) 

^In [171 section 4.1] the total derivative was denoted 5x- This is misleading as Dx = 
^(ui ,1)1 ) ■ Also Sq is the unique C-linear derivation on A. satisfying Sq{u,v) = Q and Sq o = 
Dx o Sq. 
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with if*^'" — {aun, bvn) and symmetries of similar form S = + + ■ ■ ■ with 
5*0 = 5"^'° = (cMm,(it;m)0 Here the dots may contain terms with total grading 
> 1. Certainly we have £{K°'^)S°'° = 0. The symmetry conditions with total 
grading 1 are 

£(^0,0)51,0 + £(^1,0)^0,0 (6) 

An implicit function theorem which, under certain conditions, guarantees that 
5* is a symmetry of K if the first few symmetry conditions of total grading 
0, 1, . . . are fulfilled, was given by Sanders and Wang [34l [35]. In this paper 
we restrict ourselves to solving equations (|6|. Thus we classify equations that 
admit infinitely many approximate symmetries of degree 2, which is a necessary 
condition for integrability. In the sequel we omit the adjective 'of degree 2'. 



4 The Gel'fand— Dikii transformation 

Comparing the Leibniz rule and Newton's binomial formula. 



i=0 1=0 



we see that differentiating a product is quite similar to taking the power of a 
sum. On the right hand side the index, counting the number of derivatives, 
gets interchanged with the power, while on the left hand side differentiation 
becomes multiplying with the sum of symbols. Of course, with expressions 
containing both indices and powers, one has to be more careful. The Gel'fand- 
Dikii transformation |10j provides a one to one correspondence between A^'-' 
and the space C*'-': polynomials in C[xi, ... ,Xi,yi, ... ,yj] that are symmetric 
in both the x and the y symbols. One may deduce the general rule from 

ry^ /y'^ I /-vi 1 ^2 , .3 

■jj 1 i(j o \^ u> o ^ 1 y 1 — ■ — ■ — 

U1U2V3 > — = U1U2V3, 

or consult one of the papers |T7l ETJ [24j . All usual operations from differential 
algebra translate naturally. In particular]^ 



\ G2iAaM 



5^ J, 



^We remark that only if a = fe then S may also contain terms 5^^'^^. In this paper we 
implicitly assume this does not happen. 

^As a correction to [HI Section 4.3], when (/, 9) e g'--' then / e y^'+^J and g 6 
One should think of (/, g) as representing the vector field fdu + gdv. 
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where the so called ^-functions are given by 



gl:^[a,b]{x, y) = a{x- + ■■■ + + b{y^^ + ■■■ + y^) 
-a{xi H h Xi+i + 2/1 H h 

and 

g^^^^[aM{x,y) = g{;M{y,x). (7) 

Symbolically we can solve the symmetry conditions of total grading 1 , equations 
([6|, as follows. We may write the components of the quadratic parts of S as, 
with fc = 1, 2, 

Equation ([5| has an approximate symmetry at order m with linear coefficients 
c, d iff for all z + j = 1 and k = 1,2 the right hand side of equation ^ is either 
polynomial or undefined (0/0). 



5 Nonlinear inject ivity 

In our classification we distinguish between equations whose approximate sym- 
metries necessarily have non-vanishing linear part and equations that allow 
purely nonlinear approximate symmetries. 

Definition 1 Let have total grading 0. We call nonlinear injective if 
L{K^)S — implies that S has total grading 0. And, we call an equation 
nonlinear injective if its linear part is nonlinear injective. 

With K'^ — {aun,bvn), the fc-th component of C{K^)S'^'^ , with non-zero S'*'-', 
vanishes iff Q]l''n[o.,h\ = 0. Solving the later equation with i + j = 1 yields ah = 0, 
n > 0, or n = 1, or (a — 26)(2a — 6) = 0, n = 0. In Table [l] we have displayed 
all and corresponding 5^, such that the equation {ut,Vt) = K'^ + K^, with 
arbitrary G g^, has purely nonlinear approximate symmetries 5^ G g^. Note 
that the classification is performed up to linear transformations. In particular 
we may interchange u and v. Therefore without loss of generality we set b = 1 
and classify the values of a up till inversion. 





{0,v) 


{2u,v) 


{aui,vi),a ^ 1 




{0,Vn), n> 1 




gl,0 


0-1.2 


^^2,0 ^ ^0,2 







Table 1: List of K° and such that C{K^)S^ = 0. 
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For the same choices of and the linear equation {ut,Vt) — has 
symmetries {cum,dvm) + for all m € N and c, d £ C. Indeed, every B- 
equation, that is, an equation of the form ^ with e 0^^'^, admits the 
zeroth order symmetry (2u,v). In fact, every tuple e g is a symmetry of 
{ut, Vt) = (uijVi). Or, in other words, (ui, Vi) is a symmetry of every equation. 

Only a subset of the equations {ut,vt) = K'^ + K^, with particular G 0, 
has infinitely many symmetries with non- vanishing linear part. There is a good 
reason for including such equations in the classification: their approximate sym- 
metries may correspond to approximately integrable nonlinear injective equa- 



tions. One integrable example, equation ( |37[ ), is given in section 11 On the 
other hand, nonlinear injectivity is one of the conditions in the implicit function 
theorem of Sanders and Wang, see section |3] 



6 Necessary and sufficient conditions 

In this section wc introduce convenient notation, we give necessary and sufficient 
conditions for a nonlinear injective equation to be approximately integrable, and 
we outline how we perform the classification. 
The components of equation ([s]) are 

I r^l-O I r^O.l , 7^-1.2 , \ 
bVn + + 2 + 2 + • • • / 

We denote the symbolic representation of the 6-tuple kI''^ , K'I'^, K^^'^, ^2'^, 
ii^'", K2 ^^ by . And similarly we write S'^'", . . ., 82 ^^ [> and Qn = 
5]^;°, . . . , G2-n^- ^ 6-tuple H is called proper if it consists of polynomials with the 
right sjonmetry properties, that is, if F G C2'°(g)C^'^(8)C°'2(8)C°'2(g)C^'^(8)C^'°. 
Thus, K^, S^, and ^„[a,&] are proper tuples. We will also consider s-tuples, with 
s < 6. It should be clear from the context in which space a proper s-tuple lives. 
We say that an s-tuple H = H[i] , . . . , H[g] divides an s-tuple P = P[i] , . . . , P^g] 
if I P[i] for all 1 < i < s and we write P/H = P[^ , • . . , P[s] /H[s\ ■ We 
are now able to state the following: Equation ^ is nonlinear injective and has 
an approximate symmetry of order m with linear coefficients c, d iff the 6-tuple 
= Qm[c.,d]K^ IQn[a,b] is proper. 

Let iJ, 5m[c, d] be proper s-tuples. With m{H) we denote the set of all 
TO G N such that there exists c, d G C for which H \ Gm[c,d]. And, the set of all 
proper s-tuples H with infinite m{H) will be denoted W*, or simply H when it 
is clear from the context what s is. We organize H E H by the lowest order n 
at which H divides a 5„-tuple. By Tin we denote the set of all proper tuples H 
with infinite m{H) whose smallest element is n. 

We have the following lemma. 

Lemma 2 Equation ^ is nonlinear injective and approximately integrable iff 
there is a proper 6-tuple H with m(H) infinite, such that ^„[a,6] divides K^H . 

Proof: 
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The fact that 5„[a,6] divides a proper tuple implies that equation (|9| is 
nonlinear injective. The equation is approximately integrable because for 
every m £ m{H) there are c, d such that 

^1 ^ GrnicA K'H 



is proper. 



^ Because equation ([9| is nonlinear injective, the tuple — Q„i[c,d]K^ /Qnio-M 
is well defined for all m. The integrability implies that is proper for 
infinitely many m G N and c, d £ C. This only happens when CJ„ = HP 
factorizes such that P \ and m{H) is infinite. □ 

According to Lemma [2] to classify approximately integrable nonlinear injec- 
tive equations it suffices to determine the set Ti^ of all proper 6-tuples H with 
infinite ni{H). This will be done using results from number theory, provided in 
section [Tj In section [s] we determine the proper divisors H G H} of infinitely 
many functions for possible i, k. And, in section joj we determine the 

proper divisors H G Ti.^ of infinitely many 2-tuples Q\\^^ , Qk\a "^^ere i ^ j 
if fc = 1. From those results we determine the set = [J^^z^'H^i in section 



10 For each n G N the set Tin is related to the set of approximate integrable 
equations at order n, which are not in a lower order hierarchy. 

We would like to provide an explicit but minimal list of approximate inte- 
grable equations from which one can derive all approximately integrable equa- 
tions. The following observation is useful. Let P and Q be proper tuples. From 
Lemma [i] it follows that if equation with = P, is approximately in- 
tegrable, then the same equation, but with = PQ, is also approximately 



integrable. Therefore, the classification in section 10 describes the divisors that 
have maximal degree. And the corresponding list of equations only contains 
equations with quadratic parts of minimal degree. 

From the results of sections |8]|9] it follows that 7i„ is non-empty for all n e N. 
That means there are new approximately integrable equations at every order. 



In section 10 we classify the highest degree divisors in Tin completely, that is, 
for any order. We are not able to explicitly list all corresponding equations, as 
this paper is bound to be finite. In section [TO] we do provide a complete list of 
approximately integrable equations of order n < 5. 

We explicitly provide the linear parts {cUm, dum) of all the symmetries of the 
equations in our list. This enables one to calculate any approximate symmetry 
in principle. This can be done using Maple code provided at [15J. We remark 
that if one multiplies the quadratic tuple of an equation with a proper tuple, 
the resulting equation may have more symmetries than the original one. As we 
will now illustrate it may also be in a lower hierarchy. 

From Lemma [2] we know that if _ff G 7i„ and C/„[a,6] | K^H, then equation 
([qJi is approximately integrable with approximate symmetries at (higher) order 
m € m(H). The following lemma applies. 
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Lemma 3 Suppose H G 7i„ and Qn[<i,b\ divides K^H. Then equation has 
more symmetries than the ones at order m £ m{H) iff there is a divisor Q S 
Ti-kKn of H, with m{H) smaller than and contained in m{Q), such that Qn[ci,h] 
divides K^Q. 

Proof: Given a divisor Q £ Hk of H such that tj„[a,6] | K^Q, it is clear that 
equation ^ has a symmetry at every order m g m{Q) with 

^1 ^ gm[c,d] K^Q 
Q Gn[aM' 

To see that the converse holds, let Y denote the set of orders of approximate 
symmetries, with m(H) smaller than and contained in Y. We need to prove 
that there is a Q such that Y — m{Q). Take m G Y\m{H) and write Qn = HP. 
Since Gn \ K^H we have = PR. The tuple = Q^K^ /Qn = QmR/H is 
proper. Since m ^ m{H), H does not divide Qm- There is a proper divisor Q 
of H such that Q \ Qm and H/Q divides i?, that is, Qn \ K^Q- Since Q \ H the 
set m{Q) is infinite. □ 



Remark 4 One can start with an equation that is not nonlinear injective, mul- 
tiply its quadratic tuple, and end up in the hierarchy of an nonlinear injec- 
tive equation. For example, apart from certain purely nonlinear symmetries, 
equation 1.2 has approximately symmetries with linear part (cUm, dvm) for any 
c,d € C when m is odd. By multiplying its quadratic tuple with the tuple 
[0, (/i2:i + /2yi)//, (yi+y2)/2,0, {iixi + i2yi)/i, {xi-\-X2)/2] we obtain the equa- 
tion 

f Ut \ ^ f aui + fiuiv + /2MW1 + gvvi \ 
\ Vt J \ vi -\- iiuiv -\- i2uvi -\- juui y' 

which has approximate symmetries at all orders m > for any c,d £ C, and, 
it is in the hierarchy of an equation of the form 0.3 iff fi — 12 ~ 0. In this 
paper we do not explicitly describe all symmetries of all approximately integrable 
equations that can be obtained from our list. 



7 Results from number theory 

Generally speaking, progress in classifying global classes of evolution equations 
has been going hand in hand with applying new results or techniques from 
number theory. For the classification of scalar equations the new result was 
obtained by F. Beukers, who applied sophisticated techniques from diophantine 
approximation theory . The Skolem-Mahler-Lech theorem, stated below, first 
appeared in the literature in connection with symmetries of evolution equations 
in [2]. Beukers, Sanders and Wang used a partial corollary of this theorem to 
conjecture that there are only finitely many integrable equations (|9| with = 
[0,0,1,0,0,0]. Their conjecture became a theorem in |3], where an exhaustive 
list of the integrable cases was produced using a recent algorithm of C.J. Smyth 
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ID, that solves polynomial equations f{x,y) = for roots of unity x,y. And, 
the classification of yB-equations was due to results on diophantine equations in 
roots of unity, again proved by Beukers \X2^ . 

However, as it turns out, we do not need entirely different results or tech- 
niques from number theory to globally classify two component evolution equa- 
tions, with homogeneous diagonal linear part, admitting infinitely many approx- 
imate symmetries. 

7.1 The Skolem— Mahler— Lech theorem 

A sequence {Um, fn G N} satisfies an order n linear recurrence relation if there 
exist si, . . . , s„ such that 

The general solution can be expressed in terms of a generalized power sum 

fe 

i=l 

such that the roots are distinct and non-zero, and the coefficients Ai{m) are 
polynomial in m. By definition the degree of Um is d — di, where di is 

the degree of Ai{m). It can be shown that the order of the sequence equals 
n^k + d [32 Q 

A generalized power sum vanishes identically, — for all m, precisely 
when all its coefficients vanish as polynomials in m, Ai(rn) = for all i. We 
prove this by induction on the degree. For d = the statement is plain, the 
functions h —> a'^ are linearly independent for distinct ai. Let S : f{m) — s- 
f{m + l) be the shift operator. Suppose d > 0. Then for some i we have di > 0. 
The generalized power sum Vm — {S — ai)Um has degree d— 1. By the induction 
hypothesis we have, in particular, ai{S ~ l)Ai(h) — 0. Since ^ this implies 
di — Q and hence we are done. 

Theorem 5 (Skolem— Mahler— Lech) The zero set of a linear recurrence se- 
quence {m e N : Um = 0} is the union of a finite set and finitely many complete 
arithmetic progressions. 

Note that an arithmetic progression p is complete if p = {f + gh : h E N} for 
some remainder / G No and difference g > /, 5 G N. Theorem[5]was first proved 
by Skolem for the rational numbers [37^, by Mahler for algebraic numbers [T^, 
and by Lech for arbitrary fields of characteristic zero [18 . The proofs rely on 
j>-adic analysis and consist of showing the existence of a difference g G N such 
that every partial sum, with < f < g, 

k 

Uf+gh - Y.(Mf + 9h)a{){c4f (10) 

4=1 

*In |32| one should replace equation 2.1.2 by equation 1.3 from |33| . 
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either has finitely many solutions h or vanishes identically. We refer to [271 155] , 
and references in there, for sensible sketches of a proof. 



If ( 10 1 vanishes identically the sum on the right breaks up into disjoint pieces 
/ C { 1 , . . . , m} each of which vanishes because the roots af , i G /, coincide 
and the sum of their coefficients ~^ 9h)ct{ vanishes identically as a 

function of the variable h. Since Ai{f + gh) does not vanish identically, each 
piece contains at least two terms. In particular, the following will be usefuU. 

Corollary 6 // the equation 

aia^ + a20i2 + • • • + flfcQj, = 0, 

with nonzero ai,ai G C has infinitely many solutions, the set {ai, q;2, . . . , a^} 
partitions into a number of disjoint subsets, such that each subset has at least 
two members, and the ratio of any two members of a subset is a root of unity. 

For instance, when fc = 3 the triple ai/a2, 02/(^3, cti/a^ consists of roots of 
unity. 



7.2 Diophantine equations in roots of unity 

The following theorems are of crucial importance for the classification problem 
considered in this paper. 

Theorem 7 (Beukers) Take to > 1 integer. Let /i, v be distinct roots of unity, 
both not equal to 1, such that v ^ fi^^ when m is odd. Then 

(1 - - = (1 - ^")(1 - i/)™ (11) 

implies /i™ ~ v"^ = 1. 

Theorem 8 (Beukers) Take to > 1 integer. Let /i, i' be distinct roots of unity, 
not both equal to 1, such that v ^ fi^^ when m is even. Then 

(l + z.'")(l-M)"-(l + M"')(l-i')™ (12) 
implies ^™ — v™ = — 1. 

Theorem 9 (Beukers) Take m > 1 integer. Let fi, v be roots of unity with 
fJ. ^ I. Then 

(1 + - /i)" = (1 - ^'")(1 - (13) 

implies ^™ = — i^™ = 1. 

Whereas the Skolem-Mahler-Lech theorem implies that certain ratios are 
roots of unity for the equation to have infinitely many solutions, the above 
theorems tell us precisely what the solutions are. In particular, they imply that 
the zero sets consist of complete arithmetic progressions only. 
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Theorems 7|8|9 are slightly more general than [TH Theorems 22,25], which 



were proved by F. Beukers. We won't repeat their proofs here, however we 
do indicate the difference between the two sets of Theorems, which is twofold. 
Firstly, in Theorems |7|8|9| we do not assume that /x, 7^ — 1. In certain cases 
this follows from [12, Proposition 24], in others one has to rely on the following. 

Proposition 10 (Beukers) If v is a root of unity such that 

(1 + zy™)2™~i = (1 - v)"", (14) 
then V = —1 and m is even. 

Proof: By Galois theory we may assume that v = e^'^*/". Taking n — \ does 
not give any solutions. If n = 2 then m has to be even. We will show there 
are no solutions with n > 2. When m = 1 there is no root of unity such that 
1 + 1^= 1 ~ v. Taking m = 2 it follows that n — 2. So we may assume that 
m > 2. 

Since v \\ + v™\ does not vanish and we have |1 + J^™| > sin(7r/n). Also 
we use |1 — J'l < 211 /n. This gives 

(2-)" > 11 - v\"' 11 + j/™|2"-i > sin(-)2"-i. 
n n 

Division by 2"^Tr/n yields (taking n > 2) 

(-)™-i>sin(^)f >.41, 

which implies (taking ni > 2) that ir/n > .64, or n < 5. When n = 3, |1 + v"^] 
equals 1 or 2, and |1 — I'l = \/3, whose m-th power does not equal 2™ or 2™~^. 
When n = 4, |l + ;y™| equals or a/2 or 2, and |1 — = V2, whose m-th power, 
with m > 1, does not equal or ^/22"^~^ or 2™. □ 

Secondly in the proofs wc do not in general need v ^ fJ- and v ^ l//i. And, 
we also note that in [12^ Theorem 25] it was mistakenly supposed that 7^ —1. 
This should have been /i" ^ =f1 depending on the sign in 12j equation (10)]. 



8 Homogeneous quadratic parts 

In this section we determine the proper divisors of infinitely many 1-tuples 
Gm = Gk]^'' for all possible choices of i, k. 

Due to equation ([7| we may take k — 1; equations of the form {ut,Vt) = 
{aUn, bvn + K) are related, by the linear transformation u <-> w, to equations of 
the form (ut, Vf) = {au^ + bv^). We start with the simplest case i = I. 
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8.1 Classifying approximately integrable scalar equations 

The Lie derivative of the quadratic part of a possible scalar symmetry with 
respect to the linear part K'^ = u„ of a scalar equation Ut = + + • • • is 
symbolically given by C{K^)S^ \> G]^S^ with C/-function 



Thus the case z — A; = 1 is equivalent to the scalar problem, which is easily seen 



by taking u = 0. The function is also proportional to Q]l\ '[a, a] so the results 



apply to the case a = 6 as well. 

In the classification of scalar equations [531 a different route was taken than 
the one we take. Namely, whereas we perform our classification with respect 
to the existence of infinitely many (approximate) symmetries, Sanders and 
Wang performed their classification with respect to the existence of symmetries 
(finitely many or infinitely many) . They showed in particular that there are no 
scalar equations with finitely many generalized symmetries, which confirms the 
first part of the conjecture of Fokas : 

If a scalar equation possesses at least one time-independent non- 
Lie point symmetry, then it possesses infinitely many. Similarly for 
TV-component equations one needs N symmetries. 

We note that the conjecture of Fokas does not hold inside the class of K-equations 
[16]. In their classification Sanders and Wang relied on the following 'hard to 
obtain' result from number theory, proved in [IJ. 

Theorem 11 (Beukers) Let r € C such that r{r + \){r'^ + r + 1) 7^ 0. Then 
at most one integer ra > \ exists such that Ql-^{l,r) — 0. 

In contrast, classifying the equations with respect to (approximate) integra- 
bility can be done using the following 'easy to obtain' result. Proposition [T2] is 



of course not as strong as Theorem 1 1 For obvious reasons we do not include 



the constant divisors in Tio in our lists. 

Proposition 12 The proper divisors of infinitely many Q\'^-\c,d\(l^ y) are prod- 
ucts of 

1- y E 1-12, m > 1 

2. (1 + y) e Hs, m = 1 mod 2 

3. 1 + y + y^ eH5, 171= 1,5 mod 6 

4. {l + y + y^)^ e Hr, m = 1 mod 6 

Proof: According to the Skolem-Mahler-Lcch theorem, see Corollary if 
the diophantinc equation Q^^{l,r) = has infinitely many solutions m, then 
r = 0, — 1 or r and r + 1 are both roots of unity, in which case r is a primitive 
3-rd root of unity. The orders are found by substituting the values for r. We 
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have g^,{l,0) = for all m, ^^}+2h(l'-l) = 1 + (-1)^ = ^hcn / = 1, and, 
with 1 + r + ^ 0, g}+ehi^,r) = 1 + - (1 + r)-^ = when / = 1 or / = 5. 
Fmally, by solving the simultaneous equations 0^(1, r) — drg}n{l,r) = we 
find that r is a double zero when both r and 1 + r are {m— l)-st roots of unity. 
□ 

As a particular corollary of Proposition [12] we have the following. Equation 
([9| with a = b and n — 2,3,5,7 is approximately integrable for arbitrary K^. 

8.2 i3-equations 

The case i = — 1 has been globally classified with respect to integrability in [T^]. 
This class of equations is particularly nice because any approximate symmetry 
is a symmetry. We go through the main ideas and formulate the results slightly 
different from |12| . minimizing the role of biunit coordinates. This makes the 
argument cleaner and sets the stage for the main results of this paper. As the 
case c = d is covered in the previous section, this will be excluded in what 
follows. 

Proposition 13 All proper divisors H o/t/]T^'^[c,(i](l, y) with d and m{H) 
infinite can he obtained from the following list. 

1. l + v eUi, m=\ mod 2, d 7^ 

2. (1 + y)" e Un, m>n,d=0 

3. {y - r){ry - 1) G 7^2; r ^ -1, rn > 1 

4- (y — r)^{ry — 1)^ G 7i„, r 7^ —1, n > 3 the smallest integer such that 
r"~^ = 1, m = 1 mod n — 1 

5. {y - r){yr - l){y - f){yf - 1) G Tin, r = - - 1), ^,1/ roots of 
unity such that (/i — l)(i^ — 1)(^ — — 1) 7^ 0, n > 3 the smallest 
integer such that /i" = i^" = 1, m = mod n 

6. \ + y" G Hn, m = n mod 2n, c—Q 

Unless stated otherwise, the coefficients of the linear part of the symmetries 
satisfy c/d^ (1 + r'")/(l + r)™. 

Proof: We study the zeros of the function 

g-.'^\,d]{\,r)^d{\ + r-)-c{l + rr. 

Take d 7^ 0. Then r 7^ —1 is a zero when 

c 1 + r™ 

d " (1 + r)"' ^^^^ 
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in which case 1/r is a zero as well. The point r — ~1 is a zero when m is odd, 
where it has multiplicity 1, or when d = 0, where the multiplicity is m. 

The other multiple zeros are obtained from setting the r-derivatives of the 
function to zero, see also [5]. Taking r 7^ — 1 and solving the simultaneous 
equations Gi-m^il, r) = drGi-lfi^, r) = yields r™~^ = 1, while drGi-ll^il, r) = 
d'^Gi-ll^i^, r) — yields r — —1. Therefore, all multiple zeros r 7^ — 1 are double 
zeros. We have c/d — 1/(1 + r)"'^^ and 1/r is a double zero as well. There are 
no other double zeros since the equations \r\ = \s\ and |1 + r| = |1 + sj imply 
that r = s or r = s. Let n be the lowest integer such that r"^^ = 1, so r is a 
primitive (n — l)-st root of unity. All m such that r™~^ = 1 are m = 1 mod 
n- 1. 

To classify higher degree divisors we have to find all r, s e C, with (l + r)(l + 
s){r — s){rs — 1) 7^ such that the diophantine equation 

Um{r,s) ^ G^^^[l + r"\{l + rr]{l,s) 

= (1 + r)" + ((1 + r)s)" - (1 + sr - ((1 + s)r)™ = 

has infinitely many solutions m. According to the Skolem-Mahler-Lech theo- 
rem, see Corollary |6] either rs — or one of the pairs 

1 + r (l + s)r 1 + r 1 + s 

1 + s (l + r)s r(l + s) \s(l + r) ^ ' 

consists of roots of unity. When rs = we have c — d which we exclude. 



Suppose the first pair of ([l6| consists of roots of unity, let ^— (1 + r)/(l + s) 

We may write 

M{^Ji, v) = u- 



and V = (l + l/s)/(l + l/r). We may write r = where 

M-1 



i^-l' 

and find that s = A4{1 /fi, l/v) = f. In terms of roots of unity /i, v we have 

U^n.{r, s) = [ ^^^^ly ) ((1 - m)™(1 - - (1 - - M")) • 

Note that l)(/i — v){^v — 1) 7^ because (r — s){rs — 1) 7^ 0. Hence, 

using Theorem [7] we obtain /i™ = v"^ = 1. When n is the lowest integer such 
that n"' = = 1, fi or v are primitive n-th roots of unity. And all m such that 



_ j^m _ ^ g^j.g gjyg]^ by m = mod n. Next, suppose the second pair of ( 16 1 
consists of roots of unity. By a transformation r —> 1/r we get the first pair. 
Since A4{ij,,i/)~^ = M{l/v, l//i) we get the same solutions, but with s = 1/r. 
Finally, when r, s are roots of unity Urnir,s) — can be written in terms of 
= — r, ly = — s, 

C/„(r, s) - (1 - m)™(1 + i-uD - (1 - i^ni + {~f,r). 
When m is odd Theorem |7] applies and when m is even Theorem [8] applies. □ 
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Consider the set of points 
{reC:r = z^), /x" = = 1, (a^ " ^){'^ " " - 1) 7^ 0}. (17) 

To illustrate where these points lie in the complex plane we use biunit coordi- 
nates. Suppose ipy(t' such that {ipl — \(j)\ — I and r is the unique intersection 
point of the lines ipM. and (pM. — 1. Then r = TZ{ip, (p), with 

and (j)) are called the biunit coordinates of r. Denote further 

7^(A, B) = {reC:r = 7^(a, b),ae A,b£ B, ^ b^}, 

and 

$™ = {reC:r" = l,rVl}- 
Using the algebraic relation A^(0^, = Ti{ilj,4>) one verifies that the set 



(17) is equal to {r € 'R-i^2m,^2m) ■ \r\ ^ 1}. For n = 7 the upper half of this 



set is plotted in Figure [T] 




Figure 1: The corner points r, with |r| 7^ 0, 1 satisfy Umif,''') — when rn = 
mod 7. The two circles are Irl = 1 and |r + 11 = 1. 



8.3 Quadratic terms bilinear in u-, and f-derivatives 

This section deals with the case i = 0. 

Proposition 14 If H is a proper divisor 0/ ^5'.'^[c,(i](l, y) with m{H) infinite 
then H is a product of the following polynomials. 

1. y <eHi, ■m> 0, c^O 
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2. y'^ eUn, m>n, c^O 

3. (y - r) e 7^2; r =^ 0. m > 1 

4. {y-r){y{l + r) + r) eHs, ni=l mod 2. Mien (1 + r)^" = 1, r 7^ 0, we 
also have m = mod 2n, d = 0. 

5. {y — r)^ e 7i2n; r 7^ 0, n </ie smallest integer such that (1 + 7-)2n-i _ 
TO = 1 mod 2n — \ 

6. {y — r)^(y(l + r) + r)^ G W2n+i, n > 1 the smallest integer such that 
(1 + r)2" = 1, TO = 1 mod 2n 

7. {y-r){y{l+r)+r){y-f){y{l+f)+f) G n > 3 odd, r = 

(/i — — ^ ^ 1) 7^ 0; smallest integer such that /i" = 

J/" = 1, m = n mod 2n 

- - »~) e 7Y„, n > 2 even, r = {^i - v) / {v - 1), - l){v - l)(/i - 
v){^ii' — 1) 7^ 0, n i/ie smallest integer such that fi" = i^" — 1, to = 
mod n 

9. {y~r){y{l+f)+f) G Hn, n> 2 even, r = {v - ^i) / {^i-l) , {n-v){^u-l) ^ 
0, n the smallest integer j or which /i" = i^" = —1, m = n mod 2n 

Unless stated otherwise, the coejficients of the linear part of the symmetries 
satisfy c/d = r™/((l + r)™ - 1). 

Proof: We are after the zeros of infinitely many 

c;°iM](l,y) = c-c(l + yr + dy™. 

Take c 7^ 0. Then r 7^ is a zero of precisely when 

d (1 + r)™ - 1 



c 



(18) 



When TO is odd — r/(l + r) is a zero as well. The point r = is a zero for all 
c,d,m. It has multiplicity 1, except when c — where the multiplicity is to. 
One can show that the multiple zeros r 7^ of ^".^[cjfi] are the double zeros 
{r 7^ : (1 + r)™~^ = 1}, with c/d ~ r'"^^. When r is a double zero the only 
other double zero is r = + when to is odd. 

Higher degree divisors are given by distinct non-zero r, s G C, with r+rs+s 7^ 
when to is odd, such that the diophantine equation 

U,n{r,s) = g",:l^[r"^,{l + rr-l]{l,s) 

= r™-r™(l + s)'" + s'"(l + r)'"-s" = 

has infinitely many solutions to. The cases r = — 1, s = — 1 yield the primitive 
third roots of unity, as in Proposition |12[ where c — d — 1, which we excluded. 
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Then, according to the the Skolem-Mahler-Lech theorem, at least one of the 
pairs 

r r(\ + s) s s , ^ , ^ 

1 ^ or . ,-(1 + or l + r,l + s (19) 

s s(l + r) r(l + s) r 

consists of roots of unity. Suppose the first pair consist of roots of unity. Let 
fi = r/s and jy = r(l + s)/s/(l + r). Then (^ - - l)(/i - i/) ^ 0, /zj^ ^ 1 
when m odd, r — JV{fi, v) and s — A/'(l//i, ^/v) — r with 

AA(m,^) = ^. 

u — 1 

When fj,iy — 1 and m even we have r = — r/ (1 + r). In terms of /i, ^' we get 

Um{r, s) = ((1 - /i)™(l - ^'") - (1 - ^)'"(1 - m")) , 

which implies, using Theorem [7] that /i™ — v"^ = 1. In bi-unit coordinates we 
have r S 7?.(<i>2m, <i>2m) such that |r + 1| 7^ 1 when m odd. 

Next, suppose that the second pair of (19) consists of roots of unity, ji = 
-r/s/(l + r), V = -(1 + s)r/s. We have (/x - l){v - l)(/i - i^)(M'^ ^ 1) 7^ 
when r + rs + s 7^ 0, that is, when m odd. When r + rs + s = and m even 
we get (1 + r)™ — 1, which corresponds to 6 = 0. Otherwise, r = JC{n^v) = 
(i^ — fJ-)/{fJ- — 1) and s = — r/(l + f). In terms of /i, we have 

^"^"'"^^ (mm^^T)^^) ((i-^r(i + (-M)")-(i-Mr(i + (-^r)) 

When m is odd Theorem [7] implies /i™ = i/™ = 1 , while for m even Theorem 
[8] yields ^™ = jy'" = —1. The biunit coordinate description can be found as 
follows. Solve the simultaneous equations JC{iJ.,v) = ??.(?/', </>), /C(l//x, 1/;/) = 
TI{1 / ip , 1 / (f)) to find that n = ilP/(jP,v = ifP' . For odd m we don't find new 
values for r, but for m even we get r S TZ{'^^rn \ *i'2m , *i'2m) such that |r + 1 1 ^ 1. 
Finally, suppose that the last pair of ( 19 1 consists of roots of unity. Then 
jjL = 1 + r and v ~ \ + s satisfy equation (111. According to Theorem [t] we have 
(1 + r)™ = (1 + s)™ = 1, that is, the second eigenvalue equals 0. □ 



Actually, when m is odd the two cases i = — 1, i = are related. We have 
C?°i[c,d](l,r) = giT^;'[c,d](l, -1 - r). (20) 

Indeed, at odd order m the zero r = — 1 of Qil^ translates into the zero r = 
of Q^'ln- Also the image of the unit circle jzj = 1 under /s : r ^ — 1 — r is the 
unit circle I2: + 1| = 1, relating the double zeros of the two tj-functions. The 
symmetry f2 '■ r 1/r is translated into /4 = /s o /2 o /s : — *■ ~7'/(l + f)- 
And we note that set TZ(^rn,^m) is invariant under the group of an-harmonic 
ratios, generated by /2 and /a, cf. |5U]. Using the above, for odd m one may 
obtain Propostion [14] from Proposition [13] and vise versa. 

Summarizing this section, it implies that equations with homogeneous qua- 
dratic parts are approximately integrable when n < 4. At any order n > 4 a 
finite number of new approximately integrable equations has been found. 
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9 Non- homogeneous quadratic parts 



This section deals with equations whose quadratic part is not homogeneous, that 
is, = Kl'^~^ , KI'^~-' with i ^ j when k = 1. We provide the corresponding 
sets Tin of 2-tuples. This time we do find conditions on the ratio a/b for low 
orders n < 4. 

When i = 1 the first part of the condition H G Hn, ^^[i] G Ti.^ being a 
divisor of infinitely many Gi-m^ does not give conditions on c/d, see Proposition 
In this case the are obtained from the classification of H[2] € dividing 
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infinitely many GI'.Iji'', which was obtained in the previous section. A similar 
remark can be made when (j, fc) — (0,2). Due to equation Q there are four 
cases left to consider, with fc = 1: ~ (—1, 0); and with k — 2: — (0, 1), 
(z,j) = (-l,2), = 

There are certain divisors of infinitely many fj„i[c, dj-functions for any value 
of c/d. These will be called trivial divisors. Apart from the constant divisors 
we have 

(1 + I Gi,2ni+li'^:y) V I 

(X + 1) I Gl'a^+iix, 1) X I Glfrni^, 1) 

We may take (or TJpj) to be trivial. Then e HMf iJpj (ir[i]) is one of 
the divisors of infinitely many ^-functions presented in the previous section. In 
the sequel we assume that neither nor H\2] is trivial. Also we will assume 
that cd{c -d)y^O. 

Proposition 15 We list the non-trivial divisors H of the 2-tuple Gi-l^[cjd\(l, y), 
Gl'.^[c,d]{x , 1) with m{H) infinite. Firstly suppose n is odd and P{y) divides 
Gi.l^[c,d\{lTy) with infinite m{P) whose smallest element is n, of Proposition 
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Then P{y), P{—1 — x) G Tin- Secondly, when n is even we have: 

1. (y - r)(r?/ - 1),.T + 1 e 7i2, 7- e $3, m = 2,4 mod 6 

2. {y - r)^{ry - 1)^, a; + 1 € Ji^, r € <i>3, to = 4 mod 6 

3. {y — r){ry — l),fx + f + I & Hn, r = — I)/ ^/{v — I), ji ^ I, n the 
lowest integer such that fi" = — j/" = 1, m = n mod 2n. 

The linear coefficients of the symmetries satisfy c/d — (1 + r'")/(l + r)"*. 

Proof: When the order of the equation n is odd, no new conditions on the 
linear part are obtained since the relations ^ and (20 1 imply that 



Si;l,fhd]{'^^ -1 - r) Gl:X[c,d]{r, 1). 
For even n, there should be r G C with s 7^ such that 
Si:lnhd]i'^,r) = Gl:°[c,d]{s,l) = 0, 
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or, equivalently, 



Um{r, s) = + {rsr + (1 + r)" - ((1 + r)(l + s))™ - (21) 

for infinitely many rn including n. Then, using the Skolem- Mahler-Lech theo- 
rem, we may infer that either rs(l + r)(l + s) = or at least one of the pairs 

rs s rs s 

(l + s)(l + r)' (l-t-r)(l + s)'TT7 ^ ^ 

consists of roots of unity. When r(r + 1) = we have c = c? or c? = 0, which we 
excluded. When s = — 1 we are left with the equation Um — (—1)'" + {—r)"^ + 
(1 -I- r)'" = 0. Applying the Skolem-Mahler-Lech theorem we see that both r 
and 1 + r are roots of unity and hence, that r is a primitive third root of unity. 
One verifies that U^+3k = ((-1)* + {-rY + (1 + r)*) (-l)*^ when i equals 1 
or 2. Also, if —1 is a zero of then c/d = —(—1)™. 



Suppose the first pair of (22 1 consists of roots of unity. Writing equation (21 1 
in terms of /i = 1 + s, v ~ — r, we get equation (13 1. Theorem [9] then implies 



/x™ = —J/™ = 1, which corresponds to the case a — 0, which we excluded. 
Suppose the second pair of ( p2| consists of roots of unity. Then /Lt = s/{l + 
s)/(l -I- r) and u = — rs/(l + r) are roots of unity, and we get r = K{u,v) = 
-v{^i - l)/^i/{v - 1), s = -(1 -I- f)/f, and 

/ — \ ™ 

- ( (,_;^)(;_1)J (d + - - (1 - - nuD. 

When m is even. Theorem [9] yields /i™ = —v"^ = 1 or /i = 1. But, when = 1 
we have s — —{l + r)/r and Um — 2(1 + r)™ = iffr ~ —1, which we excluded. 
Using /C(l/02, VV</'^) = ^(Va (t>) we may write r e 7^($4„J \ $2m, ^'2m)- When 
n is even and n is the lowest integer such that /z" = — i^" = 1 we have fi is 
a primitive n-th root of unity or j/ is a primitive 2n-th root of unity and all 
solutions to /i" = —v"^ = 1 are given by to = n mod 2n. 



The third pair of (22) is obtained from the second by /2 : r ^ 1/r. Under 
this transformation we have TZ{ifj,(j)) Tl{ip~^ , Hence we get the solu- 

tions r e TZ{^4m \ $2m, ^4m \ ^2m) and s = — 1 — f. Or one can express Um ~ 
in terms of — rs/(l + r)/{l + s), v — —(1 + r)/ s to find these values. Another 
way of describing the last item would be: 3. [{y — r){ry — 1), a; -I- f -I- 1] S Tin, 
r = ^{v — l)/(/i — 1), /i 7^ 1, n the lowest integer such that /x" — —v^^ — 1, 
m = n mod 2n. □ 

In the remaining cases the diophantine equation we obtain from the zeros of 
the (J-functions will be of the form 

(l + aA'")(l + 6B™) + cC™ = 0. (23) 



Lemma 16 Suppose that the diophantine equation (23), with ABC ^ 0, has 



infinitely many solutions. Then A, B , and C are roots of unity. 
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Proof: Using Corollary [6j three of the numbers 1,A,B,AB,C have a root of 
unity as a ratio and the same is true for the remaining two. Therefore at least 
one of the pairs C, A; C, B; C/A, B; C/B, A consists of roots of unity. When C 
and A are roots of unity, their powers yield a finite number of values. Moreover, 
for the infinite number of solutions we have (1 + aA"^) ^ 0. Hence, for these 
infinite number of solutions (1 + 6i3™) has only finitely many values. This only 
happens when i? is a root of unity. The other cases lead to the same result, e.g. 
when CjA and B are roots of unity we divide the equation by and find that 
^ is a root of unity. □ 

Suppose that the triple CtVt fiC^v) consist of roots of unity. Then we can 
apply the algorithm of Smyth |3] to solve the equation f{(, r])~^ = f{C^^,ri^^) 
for roots of unity. In particular, a finite number of values will be obtained. We 
denote the set of all primitive n-th roots of unity by <I>^ . 

Proposition 17 We list the non-trivial divisors H of the tuple Q^.l^l^[c,d](l,y), 
Si'mi^iy) with m{H) infinite. 

1. ?/2 + y+l,y + l e Wa, m = 1,2 mod 3, c/d= -(-1)™ 

2. {y'^ + y+\f,y + \ eUi, m=\ mod 3. c/rf = -(-l)"* 

3. {y-r^){y-f^), (y-r+ l)(y~r+ 1) e H3, r e ^'^j, m ee 1, 3, 7, 9 mod 10 

4. {y - r'^Y {y - f'^Y , {y ~ r + \Y{y ~ f +\)'^ ^Hn, r ^ m = \ mod 10 

5. {y-r){y~f), {y - r +l){y - f +1) H^, r m = l,b,7 ,11 mod 12 

6. {y~rf{y-fY, (y - r + - r + 1)^ e 7^13, r G $12, m = 1 mod 12 
The linear coefficients of the symmetries satisfy c/d = (r — l)"y(r™ — 1). 
Proof: We have ^^^^^[c, d](l, r) = ^^.'^[c, s) ^ when 

[/„(r, s) = (1 + r™)(l - (1 + s)™) + (s(l + r))™ = 0. (24) 

We want to classify all r, s e C, with rs{\ + r) ^ 0, such that equation p4| 
has infinitely many solutions. According to Lemma |16| we have s — —1, or 
r, l + s,s(l + r) consists of roots of unity. When s = — 1 we obtain that r is a 
third root of unity and Ui+sk = iff i = 1, 2. When x — r,y = 1 + s, f = s{l + r) 
consists of roots of unity, then x, y are cyclotomic points on the curve 

l + {xy- 2{x - y)){xy - 1) + {x - yf = 0, 

and can be found algorithmically. They are a; G <i>3, y e <i>g; y G ^'iq, x — y^ or 
X = if", y G ^127 X — y or X — y. The first case only happens when c = d. In the 
second case we have f — y'^ or f ^ y^ and we find C/i+iofe = iff i G {1, 3, 7, 9}. 
Note that with s = y— 1, |y| = lwe have —s/{l + s) — s. The last case gives 
/ = y4 or / = y3 and U,+i2k = iff i G {1, 5, 7, 11}. 
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The multiplicity of the zeros is obtained from Proposition 13 and 14 We 
have r S $3 is a double zero of Gi-^^ when to = 1 mod 3. When y G we 
have that both r = and r = y^ are in $g . They are double zeros of 5^^^ for 

TO = 1 mod 5. Also, we have that s = j/ — 1 and s are double zeros of Gi-m for 
TO = 1 mod 10. A similar argument shows the multiplicity in the last item. □ 



Proposition 18 We list the non-trivial divisors H of the tuple Gi-m^[c,d]{l,y), 
G2'm^{x,l) with m{H) infinite: 

1. (y — r){y — f), [x — r){x — f ) G 7^2, r G $3, to = 1, 2 mod 3 

2. {y - r)^{y - f)^, {x - r)^{x - f)^ G Hi, r e m = 1 mod 3 



3. {y - r^){y - r^), {x - r)(a; - r) G Hs., r G 



1,3,7,9 mod 10 



4- {y - r'^Y{y - f'^Y, {x - rY{x - f)^ G Hn, r G $5, to = 1 mod 10 

5. (y + r)(?; + f), (a; - r)(x - f) G 7Y4, r G $'12, TO = 1,4, 5, 7, 8, 11 mod 12 

6. {y + rY{y + f)^, {x - rY{x - f)^ e H13, r G $'12, to = 1 mod 12 
The linear coefficients of the symmetries satisfy c/d — (r + l)"'/(r™ + 1). 
Proof: We have G^.^^[c,d]{l,r) = Gl:~![c,d]{s, 1) = when 

C/,„(r, s) = (1 + r'")(l + s™) - ((1 + s){\ + r))™ = 0. (25) 

We want to classify all s G C, with rs{s + 1)(1 + r) ^ 0, such that equation 
(25) has infinitely many solutions. According to Lemma 16 the points r, s, and 
(1 + s)(l + r) are roots of unity. Hence r, s are cyclotomic points on the curve 

1 + (rs + l)(rs + 2(r + s)) + [r + sf = 0. 

Smyths algorithm yields: r, s G $3; s G $5, r = or r = s^; s G 3'']^2' J" = — s or 
r = —s. Substitution these into the equation (25 1, we obtained, by performing 
some Groebner basis calculations, the solutions to = 1,2 mod 3, to = 1,3,7,9 
mod 10, and to = 1,4,5,7,8,11 mod 12 respectively. The multiplicities are 
determined using Proposition 13 and using relation ([t]). □ 



Proposition 19 We list the non-trivial divisors H of the tuple G'i^[c^d\{l,y), 
52'm(a;,l) with m{H) infinite: 

1. J/ + 1— r, .x + 1— rG 7^2; r = 0, TO > 1 

2. {y + \ + r^){y+l + f'^), (x+ 1 - r)(a; + 1 - f) G H3, r G $io, to = 1,3,7,9 
mod 10 

3. {y+\ + r^f{y+\ + f'^f, {x + \ - rf{x + \ - ff G Uu, r G $io, m= 1 
mod 10 
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4- (y+l + r), (x + l-r) e 7^2, r e <^>[2, m=l,2, 5, 7, 10, 11 mod 12 

5. {y+l + r){y+l + f), (x + 1 - r){x + 1 ~ f) G H5, r G ^[2, = 1,5,7,11 
mod 12 

6. (2/ + 1 + rf{y + 1 + ff, (.T + 1 - rf{x + 1 - f)^ e 71:13, r e $'12; ™ = 1 
mod 12 

The linear coefficients of the symmetries satisfy c/d ^ (r™ — l)/(r — 1)™. 
Proof: Similar to the above, 5°.'^Jc, d](l, r) = S2;m[c, d](s, 1) = when 

(1 - (1 + r)™)(l - (1 + sD - (rs)" = 0. (26) 



We want to classify all r,s e C, with rs ^ 0, such that equation (26 1 has 
infinitely many solutions. If one of r, s equals —1, the other is a third root of 
unity. When r = s = — 1 we have a/b = —(—1)™, otherwise c ^ d. Suppose 
that (1 + r)(l + s) 7^ 0. According to Lemma 16 1 + r, 1 + s, and sr are roots 
of unity. Then x=l + r,y = l + s are cyclotomic points on the curve 

l + {xy+ l){xy - 2{x + y)) + (x + yf = 0. 

They are: x, ?; € $g; y G $io, x = -y^ or a; = - f-; y e $12, a; = or x = -y. 
The first are zeros only when c = d and the others yield the results. □ 



10 Global classification of maximal degree divi- 
sors 



Combining the results obtained in Propositions 12 13 14 15 17] [18] [19] we de- 
termine the set of all highest degree proper 6-tuples with infinite m(H). This is 
equivalent to a global classification of approximately integrable two-component 
equations with a diagonal linear part, see section [6] Highest degree tuples are 
formed as follows. With H E 7i„ and F £ Hk we have HF G Hi, where I is the 
smallest number in m{HF) = m{H) n m{F). 

Clearly if H divides Gm[ci d\ then H divides Gm[c/d, 1]. Also, if 

H = [0{xi,X2),P{xi,yi),Q{yi,y2),R{yi,y2),S{xi,yi),T{xi,X2)] 

divides Qm[c, 1] then, according to equation ([t]), the function Qm[^/c, 1] admits 
the proper divisor 

H'^ [R{xi,X2),S{yi,xi),T{yi,y2),0{yi,y2),P{yi,xi),Q{xi,X2)]- 

Thus we scale d in Gm[c, d] to 1, and perform the classification up to inversion 
of c. 

We also include tuples with zero components in the list. They correspond to 
equations that are not nonlinear injective, see section [U For each K° in Table 
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[T]we have determined the highest degree r-tuple H, with m{H) infinite, which 
divides the r-tuple consisting of the non-zero components of its C/„-tuple. The 
quadratic tuple has r non-zero components = Gn/H, unless complemen- 
tary components of the fj^-tuple vanish at infinitely m e m{H). 

First we deal with Ti^^g. Here we translate our symbolic results into dif- 
ferential language. For any H G Tin of highest degree, which divides t/„[a, 1], 
we determine from — [e, f, g,h,i, j]Qn[a,l]/ H . In principal, the tuple 
[e, f,g,h,i,j] may consist of proper polynomials, it is a common factor of ii'^ 
and S^. However, when writing down the differential equation, the e, /, g, h, i,j 
will appear in it as constants, and any other constants will be absorbed by 
them. This organizes the quadratic part of the equations. At the same time 
it may remind the reader of the fact that the quadratic tuple of the equations 
can be multiplied by arbitrary proper tuples. We give each maximal degree 
tuple H £ TLn two indices, H = Hn.h, where n is the order, and ft, is a counter. 
We label the corresponding approximately integrable equation by n.h. And, we 
express the linear coefficients c/d of the approximate symmetries in terms of 
integer sequences, or in its power sum solution if that displays well. 

Secondly we give a general description of Ti-nyb- Using this result one can, 
in principle, write down the corresponding approximate integrable systems at 
any particular order. Some Maple code has been provided at |15j . 



0.1 



10.0 Zeroth order 

According to Table [T] there are two special values of a related to equations that 
are not nonlinear injective. We have the following 6-tuplcs in Tio- At a = 

i/o.i = [0,1,1,1,0*,!]. 

Both the first and fifth component of the zeroth order t/-tuple vanishes. How- 
ever, there exist higher order CJ-tuple with zero first component, but no higher 
order fj-tuples with zero fifth component. This is denoted by the *, which indi- 
cates that in the approximate integrable equation the term K2'^ vanishes. The 
equation 

ut \ _ ( ew^ -I- fuv + gv^ 
vt ) \ V + hv"^ + ju^ 

has approximate symmetries at order m = 1, for all c/d £ C, and at any order 
m, with c = 0. 

At a = 2 we have 

ffo.2 = [1,1,0*, 1,1,1]. 

The equation 

Ut \ _ f 2u + eu^ + fuv 
Vt J \ V + hv'^ + iuv + ju^ 

has approximate symmetries at any order m G N, for all c,d £ 
And, at generic values of a we have 

i/0.3 = [1,1,1,1,1,1]. 



0.2 
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The equation, with a{2a — 1) (a — 2) 7^ 0, 

Ut \ _ f au + eu^ + fuv + gv^ 
Vt I \ V + hv^ + iuv + ju^ 

has approximate symmetries at any order m € N, for all c,d 



10.1 First order 

We have the following 6-tuples in T-Ci- At a = 1 

Hj^i = [0,0*,0*,0*,0*,0*]. 

The equation 



0.3 



1.1 



has approximate symmetries at any order m G N, for c = 0. Of course, when 
e = any S" G g is a symmetry of this equation. 
At generic values of a we have 

H1.2 = [0,yi,yi + y2,0*,xi,xi + X2]. 

The equation, with a 7^ 1, 



Ut \ f aui + ev? + fuv + gv^ 



Vt I \ vi + iuv + ju^ 



1.2 



has approximate symmetries at odd orders m = 1 mod 2 for c = 0. Again, 
when e = there are more approximate symmetries, namely at odd orders 
m = 1 mod 2 for all c, d G C, see Remark |4j 

10.2 Second order 

At second order there is the tuple 

H2 = [xiX2,yiA,yiy2,xi,l\ 

which divides all higher t/m[c, d] for all c,(i G C. The maximal degree divisors 
H2.i € 'H2 are H2.i = H2T2.i with 



T2.1 = 


[0, yi,yf + yl, 1, xi + 2yi, {xi + X2f] 




T2.2 = 


2/1^2, ^,yi,XiX2] 




T2.3 — 


[l,xi + yi,yl + yiy2 + yl, l,xi + yi,xl + X1X2 ^ 


^xl] 


T2A — 


[1, 1, {ry2 - yi){ryi - ^2), 1, 1, 1] 




T2.5 = 


[l,rxi - yi, 1, 1, 1, 1] 




T2.6 = 


[1, 1, (yi + 2/2 - ty2)(2yi + y2 + y2i), l, + yi 


+ 2/1'-, 1] 


T2.7 = 


[1, 2yi - ixi + 2xi + xij, 1, 1, 2x1 - yil - i^yi + 


2yi,l], 
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where = —1,7^ = 3. 
The equation 



lit \ / eu^ + fuv + \ 2 1 

vi J \ V2 + hv^ + iuv + ju"^ J 

has approximate symmetries at orders m = 2 mod 4 with c = 0. The equation 

?ii \ ^ / U2 + eu^ + fuv + gv'^ , ^ 2 

Vi J \ V2 + hv^ + iuv + ju^ 

has approximate symmetries at all orders m > 1 with c = d=l. The equation 

Uf \ ^ / -U2 + eu^ + fuv + gv"^ \ ^ 3 

Vt ) \ ^^2 + hv"^ + iuv + ju"^ J 

has approximate symmetries at order m = 1,2 mod 3, with c/d = —(—1)'". 
The equation, with (1 + r'^){l + r) ^ 0, 

ut \ ^ f ji^U2 + eu^ + f{ruvx - (1 + r^)uxv) + gv^ \ 

'^t J y t;2 + /i^^^ + i(r-uit; - (1 + r)^uui) + j(2ra2W+ (1 + r)^uf) ) 

2.4 

has approximate symmetries at all orders m > 1 with c/d= (1 + r™)/(l + r)™. 
The equation, with r(r + 2) ^ 0, 

ut ^ / 2+f '"2 + eu^ + /wu + g(rvl - 2vv2) \ ^ 5 

vt J \ V2 + hv'^ + i{uiv+{2 + r)uvi)+j{2uu2 + {2 + r)ul) J 

has approximate symmetries at all orders m > 1 with c/d = r"*/((l + r)™ — 1). 
The equation, with — —1, 



ut \ _ f i-l + 2i)u2 + eu'^ + f{5uiV+{3 + L)uvi)+gv'^ , 



( y \ V2 + hv + iuv + j(4uw2 + (1 + ) 



has approximate symmetries, with c/d = — 1+(— l)^'" '^'>/^2"^/^l, at order m = 2 
mod 4. The equation, with 7^ = 3, = — 1, 



\ _ f t(2 + 7)w2 + + fuv + g{4vv2 + (2 + 7 - i)vl) , ^ 7 



+ fuv + g{4:VV2 + (2 + 7 - i^-"^^ 
Vt J \ V2 + hv'^ + iuv + j{Auu2 + {2 + i- j)ui) 

has approximate symmetries at orders m = q mod 12, with q G {1, 2, 5, 7, 10, 11}. 
Define integers Pk by Pi = 1, P2 = 2, and 

p _ / -Pfc-i + Pk-:i k = 1 mod 3 , , 

\Pk-i + Pk-2 k = (),2 mod 3 • ^^'^ 

When (7 = 2 or (7 = 10 the coefficients of the linear part of the approximate 
symmetries of equation 2.7 arc given by 

C ' -N(m-<;)/12 



d 
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or else by 



/_-,^,(m-9)/12C ^ f P(3m-5)/2 +^'(3m-7)/27 g-6±5 , - 

'^^ ' d 1 P(3™+l)/2 + P(3m-l)/27 g = 6±l • ^^^^ 

The sequence {Pk} is quite interesting in itself, see [38l Sequence A140827]. It 
satisfies the 6-th order recurrence Pk = 4Pfc_3 — Pk-e- And thus it consists of 
the three subsequences [351 Sequences A001075, A001353, A001835]. The first 
two of these subsequences are the denominators and numerators of convergents 
to V3, we have P|„_i - 3P|„_2 = 1. 



10.3 Third order 

At third order the product 

-^3 = H2[xi + X2, l,yi + y2,v\ + yi, l,a;i + x-^ 

divides tJm[c, when m odd, for all c,d €z C. The maximal degree divisors 
H3,i e 7^3 are Hs.i = H^T^.i where 

P3.1 = [0,yl,yl-yiy2 + yl,i,xl + 3xiyi + 3yf,{xi + x2)'^] 
P3.2 = [i,xi{xi + yi),yiy2,l,yi{xi+yi),xiX2] 

T3.3 = [^A,{yi-ry2){yir-y2)A,{rxi + {l + r)yi){xi + {l + r)yi),l] 
T3.4 = [l,xl + xiyi + {2-(p)yf,yl + (pyiy2 + yl,l,xl + {2-(l)){xiyi+yl), 

xl + {1 ~ (t))xiX2 + xl], 

in which (f> denotes the golden ratio or its conjugate, that is, (j>{(j)~ 1) = 1. Note 
that both H2A and iJ2.5l-i-r divide i?3.3. 
The equation 

( '^t \ ^ f eu^ + fuv + gv"^ \ ^ ^ 

\ Vt J \ ^3 + hv^ + iuv + jv? J 

has approximate symmetries at orders m = 3 mod 6 with c = 0. The equation 

( \ ^ ( ""3 + eu^ + fuv + gv'^ \ ^ ^ 

\ Vt ) \ V3 + hv"^ + iuv + ju^ J 

has approximate symmetries at odd orders with c = d = I. The equation, with 

f "0 = 1 .f{fuv2 - (1 - r + r'^){u2V + uivi) 

\ Vt ) \ V3 + hv"^ + iuv + j {2ruu2 + {1 + r + r^)ul) 

has approximate symmetries at odd orders m with c/d — (1 + r™)/(l + r) 
The equation, with 0(0 — 1) = 1, 

/ Ut \ _ f -{2 + 3(t>)u3 + eu'^ + fuv + gv^ \ , 
\ Vt J \ V3 + hv"^ + iuv + ju^ J 
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has approximate symmetries at order m — q mod 10, q G {1,3,7,9}, with 

where the F^, are the Fibonacci nmiibers Fq = 0, Fi = 1, = F^-i + 7^-2 [551 
Sequence A000045]. 



10.4 Fourth order 

At order four the maximal degree divisors i G 7^4 are H4,i — H2Ti.i where 





= [0, yl vl + yi 1, {xi + 2yi){x\ + 2xiyi + 2yl), {x^ ~ 


K^2)1 


T4.2 


[l,xi + 2/1, (yi + 2/12/2 + l,a;i + yi, (x^ + xia;2 


+^^)^] 


T4.3 


= [1, 1, {yl + 2j/iy2 + 2yl){2yl + 2yij/2 + yl), 1, 1, 1] 




TiA 


= [1,(22/1 + (3 -7^x1)^1, 1,1,1] 




T4.5 


= [1,2;2 + 2x12/1+22/^1, 1,1,1] 




T4.6 


= [l,2/?-2iXi2/i-(l + i)a;?,l,l,l,l] 




T4.7 


= [1, 2/? + 2(1 - ./3)xi2/i - 2i/3x?, 1,1,1, 1] 




T4.8 


= [1, 1, Ayl + (6 - 2t + 2t/3)2/i2/2 + 42/^, 1, 2xi + (1 - i 


+ t/3)yi,l] 


T4.9 


= [1, 1, 2yl + {l-L~Lp- 2/3)2/12/2 + 2j/2, 1, 2xi + (1 - 


f i + /3)yi,i] 


74.10 


= [1, 1, Ayl + (6 - 4t - 3t/3 - /3)2/i2/2 + 42/|, 1, (1 + t - 


t/3)xi + 2/1, 1] 


74.11 


= [1, 1, yl + yiy27 + 2/2> - xiX2i + xl], 




with 


= -l,/32 = 2,7^ = 3. 





The equation 



Wt y I U4 + h{Avv2 + ivl) + + ju^ 



has approximate symmetries at orders m = A mod 8 with c = 0. The equation 

4.2 



Ut \ _ / -U4 + e{Au2U + 3uf ) + f{uv2 + uivi + 2u2v) + gv^ 
Vt J \ V4 + h{Av2V + Svl) + i{2uv2 + uivi + U2v) + ju^ 

has approximate symmetries at order m = 1 mod 3 with c/d = —(—1)™. The 
equation 



Ut 
Vt 



/ — 3m4 + e(4u2U + 3it^) + /(6U3U + 9u2Wi + 6miU2 \ 

+2MW3) + gv"^ 
V4 + h{Av2V + 3vl) + i{2uzv + 2u2Vi + 3uiV2 + 2UV3) 

\ +j(4u4U + 4uiU3 + 3^2) / 



4.3 



has approximate symmetries at order m = mod 4 with c/d — l+(— 1)™/^2™/^. 
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The equation, with + C + 1 = 0, 



Ut 
Vt 



( 3(l + 2C)u4 + e(4uw2 + 3u?) + /(6uit;+(l-4C)MVi) ^ 
+g{UviV + 3(4 - Q){Av3Vt + 2,vl)) 



Vi + h{Av2V + 3wf ) + «(7m3W + (2 + 3C)(2u2'yi + 2mz;3 



V 



+3uii;2)) + j(14u4U + (2 + 3C)(4u3Wi + 3m^)) / 
has approximate symmetries at order m=\ mod 3, with 



(_3)(m-i)/2^ TO = 1 mode, 

-(l + 2C)(-3)(™"2)/2, TO = 4 mode. 



The equation, with e = ±1, 

/ e/5ii4 + e(47iU2 + 3ui) + /(2miw + (1 - 2e)Mi;i) \ 

+g{Av^Vi + Svj + (1 - 5e)viv) 
V4 + h{4:V2V + 3vl) + i{10{2u2Vi + 3uiW2 + 2uv3) 

V +{5-e)u3v)+j{5{4u3Ui+3ul) + {5-€)u4u) ] 
has approximate symmetries at order to = mod 4, with 



Vt 



1 - (-4)"/4 ■ 



The equation, with 
/ 



Ut 
Vt 



V 



t/3M4 + e(4MW2 + 3ul) + f{2uiv + (1 - 2l)uvi) \ 

+5(10w4 + (1 - 3i)(4wii;3 + 3vl)) 
V4 + h{4:V2V + 3vf) + z(5u3W + 3(3 + l){2u2Vi + 3uiV2 

+2UV3)) + j{10uiv + 3(3 + l){3uI + Ausui)) ) 



has approximate symmetries at orders m = 4 + fc8, fc G N, with 



{-\f2 



1' 



where the integers Ai are the NSW numbers defined by ^0 = ^1 ,^1 
= eAi_i - An-2 EH Sequence A002315]. The equation, with c? = -2 



/ aM4 + e(4itU2 + 3iiJ) + /(2uiv + (1 — a)uui) \ 



Ut 
Vt 



+g{2{AviVi + 3vl) + (2 + Q!)w4z;) 



V 



W4 + h{Av2V + 3'(;J) + i((l — a)u3t; + 2(2u2Wi + 3uiU2 

+2uv'i)) + j{3u\ + Au^ui + (1 — aju^v) J 

has approximate symmetries at orders m = 4 + fcS, fc G N, with 

c _ {-lf2^^a 



fe+i 
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4.8 



where the integers Bi are defined by Bo = —1 ,Bi = 1, Bi = 34i?i_i — B„-2 
[551 Sequence A046176]. The equations, with = -l,/?^ = 2, 

/ {-1 + 2l{3 + (3))u4 + e{Auu2 + 3ul) + f{6uv3 \ 

+(9 - i + 2/3)(2uiW2 + 3u2i'i + 2u3i>o)) 

ut \ +g(12w2 + (9 + /3 + Ow?) 

Vt J V4 + h(Av2V + Svl) + i{6u2V + {3 + 5l — 3l(3 — 4:P)uiVi 

+2(1 + 3t - l(3)uv2) + j{12uu4 

\ +{3 + L-2f3){3ul + Au3Ui)) ) 

( (-17+12/3 + 2t(3-2/3))u4 + e(4MU2 + 3u?) + /(10ut;3 \ 
+ (11 - 3t - 6/3 - 2t/3)(2uii;2 + 3M2t'i + ZugWo)) 
+.9(20ww2 + (17 - 2/3 - I + t/3)v2) 
U4 + h(j^V2V + 3ui) + z(10u2W + (13 + t + 7/3 + 4/3t)ui-(;i 

+2(7 - L + 5/3)mw2) + i(20uu4 

V +(9 + 3i + 2(3 + i)/3)(3u2+4u3Mi)) / 
/ (17- 12/3 + 8i(3-2/3))M4 + e(4uu2 + 3M?) + /(12uw3 \ 

-(3/3 + 4t/3 - 20 + 6t)(2uiW2 + 3u2fi + 2m3Wo)) 
+5(24w2 + (22-3/3+i/3)w2) 
W4 + ft(4w2'u + 3v\) + ^(2^2 + (1 + 2t(/3 - l))uiwi 
+(i - 5 + (4 - L)(i)u2v) + j{2Auu4 

V +(4 + 6i+(3 + 4i)/3)(3u2+4u3Ui)) / 
have approximate symmetries at orders m = 4 + fcS, fc G N, with 



4.9 



4.10 



c 

d 
c 

d 
c 



= -l + .(-l)'=22'=+i(C„+i + C,„/3), 

= (-1 + /,(-l)'=22'=+l(a„+l + anP)){C2rn+l " ^2^/3), 



= (1 + t(-l)'=26'=+3(C„,+ i + C,„/3))(C2™+1 - C2™/3), 



respectively, where Co = 0, Ci = 1, and C2n = C'2„-i + C2n-2, ^2,1+1 = 
2C2ra — C2,i-i. These integers, cf. Sequence A002965], are the denominators 
and numerators of convergents to \/2, we have 

Cl - 2Cl_i = ±1, when n = 1 ± 1 mod 4. 

The equation, with = 3 



Ut 
Vt 



f {7 + 47)7x4 + e(4uii2 + 3ul) + /(2uii;2 + 3w2i;i + 2ii3W ^^ 
+(27 - 3)mi;3) + g{6vv2 + (6 + 7)1;^) 
V4 + h{Av2V + 3u^) + i((3 + 27)u3W — 2wu3 — 3witi2 

\ -2u2Wi) + i(6-uu2 + (6 - 7)uJ) J 



4.11 
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has approximate symmetries at order m = q mod 12, q G {1,4,5,7,8,11}. 
When g = 6±2 the coefficients of the linear part of the approximate symmetries 
are given by 



where the integers Pj, are defined by the recursive formula ( 27 1 . When q is odd, 
c/d is given by equation (28 1. 



10.5 Fifth order 

At order five the maximal degree divisors -ffs.i G 7^5 are H^ ^ — H^T^j where 

H5 = H-ilxj + X1X2 + xj, 1, l,yl + yiy2 + yl, 1,1] 

and 

7^5.1 = [0,y^yf- 2/22/? + 2/?y|-y22/i+ 2/1, 1,a:i + 5x?2/i + 10a;2y2 + i0a;iy3 

+5y4,(xi+X2)4] 

T5.2 = [l,2ylxl + 2xlyi+xt + ylxi,yfy2 + yiyl + ylyiA,^ylxi + 2ylxl 

+xlyi + yf, x\x2 + x\x\ + Xix\] 
T5.3 = [I,l,ry2-(l + r2)yiy2 + r-y^,l,™2 + (l + r)2(xiyi + y?),l] 
T5.4 = [l,l,{yl + ylf,l,{xl + 2x^y^ + 2ylf,l] 

T5.5 = [l,l,2/4 + 2/2' + (2 + 0)(2/?y2 + 2/iyi) + (4 + 0)2/?yia,a:t + (2-0)(y4 

+xlyi + 2ylxi + 2ylxl),l] 

75.6 = [l,a;? + xiyi + (2-7)yi,yi +yiy27 + y2,l,2;i + (2-7)(xiyi+2/i), 

- a;ia;27 + Xj], 

with 7^ = 3 and - 1) = 1. Note that 7/5.3 = ^^sTs.a and, both H2.7 and 
iJ4.9 divide H^,^. 
The equation 

/ Ut \ _ / 67/2 + _|_ ^^2 \ ^ ^ 

\ "^t J \ 'I's + /it^^ + *™ + ju^ J 

has approximate symmetries at orders m = 5, 25 mod 30 with c = 0. The 
equation 

( ^ = U5 + eu^ + fuv + gv'^ \ ^ ^ 

has approximate symmetries at orders m = 1 , 5 mod 6 with c = d = \. The 
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equation 



/ 1 



ut 



tU5 + eu"^ + J{{r^ - + - r + l){uiV:>, + 2u2W2 



V 



{l + rf 

+2u^Vi + U4v) — r{r'^ + r + l)uv4) + g{{l + r^)u^ 

+2(r2 + r + 1)W2) 
V5 + hv'^ + + r)'^{uv2 + uivi) + (r^ + r + l)u2v) 

+j(2r(r^ + r+ l)uu4 + 2{r'^ + 3r^ + 5r^ + 3r + 1)uiM3 

+ (r4 + 5r3 + 7r2 + 5r + l)u2) / 



5.3 



has approximate symmetries, with c/d~ (1 + r"')/(l + r)™, at orders m= 1,5 
mod 6. The equation 



Ut 
Vt 



-|u5 + ev? + f{u4V + 2U3V1 + 2u2V2 + uivs + UV4) + gv'^ 
U5 + hv'^ + iuv + j(2uu4 + 6U1U3 + 5^2) 



5.4 



has approximate symmetries, with c/d ~ (— 1)*^™ i)/42(i ni)/2^ j^^^ orders m 
1, 5 mod 12. The equation, with </)(<?!) — 1) = 1, 



Ut 
Vt 



-(4 + 50)w5 + eii^ + /((4 - 0)uu4 + 11(miU3 + 2u2W2 

+2u3i;i + U4V)) + gv^ 
f5 + hv'^ + iuv + j{2(j)uu4 — 2uiU3 + (20 — 1)^2) 



5.5 



has approximate symmetries at orders to = 5,25 mod 30, with c/d = — (1 + 
Fm-i+F^cj)), where F„ denotes the TO-th Fibonacci number. And, the equation, 
with 7^ = 3, 



Ut 
Vt 



I (26 + 157)7/5 + eu^ + f{u2V + uivi + (7 - 1)^2) \ 

+.g(4zw2 + {3 + "/)vl) 

V5 + hv'^ + i{uiVi + UV2 - (1 + l)u2v) 



V 



+j(4uu2 + (3 - -f)ul 



5.6 



has approximate symmetries at orders to = 1,5 mod 6, with c/d given by 
equation ( 28 1 . 



10.6 Higher order 

Define 

Hj = Hz[xl + X1X2 + xl, 1, 1,2/1 + 2/12/2 + yl, 1, 1], 

and, for convenience, i?2n = H2, i?2n+i = Hk, where n = k mod 6 and k S 
{3, 5, 7}. We first list the divisors that have similar structure at infinitely many 
higher orders. After that we consider the exceptional cases. 
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The Ccise a = 

i. The tuple Qn[0, 1] has the following divisor in 

[0, y^,y^ + vl, X, - (xi + yi)", {xr + x^T] . (30) 
where X is the fourth component of It divides 5m [0, 1] when 

{n mod 2n, n even, or n = 3 mod 6, 
n mod 6n, n = 1 mod 6, 
n, 5n mod 6n, n = 5 mod 6. 

The Ccise n is even 

Let n > 2 be even. We have the following divisors in H^. 

ii. For any primitive (n — l)-st root of unity r 

H2[l, 1, (2/1 - ry2)'(ryi - ys)', 1, 1, 1] (31) 
divides ^^[(1 + t")^""*, 1] when m = 1 mod n — 1. 

iii. Let one of /i, be an n-th root of unity, and the other a primitive n-th 
root of unity, such that (/U — l)(i' — 1)(/^ — — I) ^ 0. Then, with 
r = i/(M -!)/(!/ -1) 

-f^2[l, 1, (yi - r-y2)(ryi - y2)(yi - ry2)(ryi - 2/2), 1, 1, 1] 

divides ^^[1 + (1 + f)"^] when m = mod n. 

iv. For any primitive (n — l)-st root of unity r 

if2[l,(yi-(r-l)a;i)2,l,l,l,l] (32) 

divides Qm[{'>' — 1)™~^, 1] when m = 1 mod n — 1. 

V. Let one of /i, u be a primitive r?,-th root of unity, and the other a n-th 
root of imity, such that (/U — l)(i' — l)(/i — J^)(m^ ~ 1) 0- Then, with 
r = zy)/(z/- 1) 

-f^2[l,(yi -rxi){yi -fa;i),l,l,l,l] 
divides Qmir"^, (1 + — 1] when m = mod n. 

vi. Let one of n,uhe a, primitive 2n-th root of unity, and the other a primitive 

2n-th root of unity which is not a n-th root of unity, such that {n—v){^v — 
1) ^ 0. Then, with r ^ {v ~ fi)/{fj, - 1) 

ff2[l, {yi - rxi)(yi(l + r) + rxi), 1, 1, 1, 1] 

divides Gm[r"^, (1 + J")™ ^ 1] when m = n mod 2n. 

vii. Let either ^ he a, primitive n-th root of unity and be a 2n-th root of 
unity which is not a n-th root of unity, or let /i 7^ 1 be a n-th root of unity 
and V a primitive 2n-th root of unity. Then, with r = —1^(11— 1) l(y— 1)//U 

1, {yi - ryi){ry2 - yi), 1,1, rxi + (1 -hf)yi,l] 

divides Gmi^ -|- r™, (1 -|- r)™] when m = n mod 2n. 
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The case n is odd 

Let n > 3 be odd. We have the following divisors in 7i^. 

viii. For any primitive (n — l)-st root of unity r define 

Qn{r) = [I,l,{yi-ry2)^{ryi-y2y,l,{xi + {l + r)yif{rxi + {l + r)yif,l]. 

The tuple Hr,Qn{r) divides + r)i-", 1] when 

f 1 mod n — 1, n = 1, 3 mod 6, 

~1 1, n mod 3(n — 1), n = 5 mod 6. 



ix. Let one of /x, be a n-th root of unity, and the other a primitive n-th 
root of unity, such that (/i — l)(y — l)(/i — v){lJ-v — 1) 7^ 0. Then, with 
r = iy{fi — — 1) we define 

Wn{r) = [1, 1, (j/i - ry2){ryi - y2){yi ~ fy2){fyi - 1/2), 1, {xi + (1 + f)yi) 
{fxi + (1 + f)yi){xi + (1 + f)yi){fxi + (1 + f)yi), 1]. 

The tuple H^Wnir) divides ^„i[l + r™, (1 + r)"] when 

n mod 6n, n = 1 mod 6, 
m = i n mod 2n, n = 3 mod 6, 
n, 5n mod 6n, n = 5 mod 6. 



We also get new highest degree divisors in 7i^, with n odd, from primitive 
p-th roots of unity with p < n. This happens when n = I mod p (or n = mod 
p) and n = k mod 6 with fc G {3, 5, 7}, such that there is no 1 < g < n with 
g = 1 mod p (or g = mod p), and q = I mod 6 with / € {3, 5, 7} and I > k. For 
example, the tuple T5 4 divides a at m = 1 mod 4. The tuple -ff5.4 = H^T^ 4 
divides a CJ„i at m = 1, 5 mod 12, see equation 5.4. We have 13 = 1 mod 4 and 
13 = 7 mod 6. There is no 1 < g < 13 such that q = 1 mod 4 and q = I mod 6, 
with Z G 3, 5, 7 and I > 7. Thus, the tuple HjT^,/^ G 7ii3 divides a at m = 1 
mod 12. We have 9=1 mod 4 and 9 = 3 mod 6, but there is g = 5 such that 
g = 1 mod 4 and g = 5 mod 6, and 5 > 3. Therefore the tuple H^^T^ 4^ is not in 
Tig. It is in 7^5 but it does not have maximal degree. 

Also, we may have p odd. If the tuple (31) divides Qm[{^ + r)^"™, 1] when 
TO = 1 mod p, then, according to Proposition 15 Q^(^^ divides ^2j9+i[(l+''')^^, 1]- 



This always give us a new highest degree tuple in Ti.2p+i- On the other hand, 
from case iii one can conclude, using Propositions [14] and [15] and inverting 
a — > 1/a, that Qp+i{—r) divides Q2p+i[{l + r)'^^ , 1]. This we knew already, since 
when r is a primitive p-th root of unity, with p odd, then — r is a 2j>-th root of 
unity, cf. case viii. 

In general one can show the following. 
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• Suppose n = mod p and n = k mod 6 with k G {3, 5, 7}, such that there 
is no (? < n with q = mod p, and q = I mod 6 with I S {3, 5, 7} and 
I > k. Then p is odd. And, 



P, 

p,5p 



p = 1, 3 mod 6, 
p = 5 mod 6. 



(33) 



Suppose n = I mod p and n — k mod 6 with A; G {3, 5, 7}, such that there 
is no q < n with q = I mod p, and q = I mod 6 with I G {3, 5, 7} and 
I > k. Then 



P- 
P- 
P- 
P- 
P- 
P- 



1 

l,2p 
l,2p 
l,2p 
l,3p 
l,2p 



l,4p 
l,3p 
1 
1 

1, 6p 



l,6p 
1 



1 



p = mod 6, 
p = 1 mod 6, 
p = 2 mod 6, 
p = 3 mod 6, 
p = 4 mod 6, 
p = 5 mod 6. 



(34) 



We can now describe all highest degree divisors at odd order ?i, involving the 
tuples Q, W . For odd n there are the cases viii and ix described above. They 
corresponds to the cases n = p + 1 in ( 34 1 , and n = p in ( 33 1 , respectively. 
Furthermore we have: 

X. If rt = 5p and p = 5 mod 6, then iJ„Wp(r) divides Qm[^ + J'™, (1 + t)"^] 
when m = bp mod 6p. 



xi. If 71 = 1,3 mod 6 and equation (34) holds for certain p, then HnQp{r) 
divides t/,„[l + r™, (1 + r)'"] when m = 1 mod n — 1. 

xii. Let n = 5 mod 6. If n = 2p+l, p = 5 mod 6 then HnQp{r) divides + 
r™, (1 + r)™] when m = l,n mod 3(n — 1). And if n = 2p + 1, p = 2 mod 
6, or n = 4p+ 1, p = 1 mod 6, then HnQp{r) divides gj[l + r™, (1 + r)™] 
when m = 1, n mod 3(n — 1) /2. 

Exceptional highest degree divisors, n >7 

xiii. The tuple Grll, 1] G Tij divides Grn[c, d] with m = 1 mod 6 and c = d = 1. 

xiv. The tuple iJrTg.g G divides ^,„[1 + r™, (1 + r)"] with m = 1 mod 6. 

XV. The tuple HjT^ 4 G Tiy divides ^„i[c, d] with m = 1,7, 13, 19 mod 30 and 
c, d given by equation |29[ 

xvi. The tuple H-jT^ q G Ti.j divides C/,„[c, d] with m = 1 mod 6 and c,d given 
by equation |28[ 



xvii. Define 



Z - [1, {xl + xm + (1 + cl))ylf, {yl + (1 - (/.)yiy2 + yjf, 1, 
(a;2 + (1 + 0)(a;i2/i + yj))^ , {xj + (j)XiX2 + xlf]. 



34 



Then H^Z G Hu divides t/m[c, cZ] with m = 1,11 mod 30 and c/d given 
by equation (29 1. And HjZ e Hsi divides t/„i[c, d] with m = 1 mod 30. 



xviii. The tuple 



Hr[l, {xj + xivi + (2 - j)ylf, (yf + 7yiy2 + 1, 
{xl + (2 - j){xiyi + yj)f, {x\ - 70:1X2 + x\f\ e H13 



divides Gm[c, d] with m = 1 mod 12 and c/d given by equation (28 1. 



11 Concluding remarks 

In the formal symmetry approach 125], as well as in the computer-assisted 
schemes HO] , not knowing the ratios of eigenvalues strongly complicates the 
classification of integrable equations. There the ratios are obtained, if possible 
at all, at the very last stage of the calculations. We hope that the a priori knowl- 
edge provided in this article will be an impetus to complete the classification. 

Usually, in classification programs one considers homogeneous equations. A 
2-component equation {ut,vt) ^ K is homogeneous of weighting A if .ftT is an 
eigenvector oi C{ax + \i(Tu-\-^2'^v), where Gx — {xui^xvi) counts the number of 
derivatives. We have compactly provided a list of nonhomogeneous equations. 
A complete list of homogeneous equations can be obtained from our list by mul- 
tiplying the symbolic quadratic parts with appropriate tuples of polynomials. 
And Lemma |3] can be used to determine all symmetries of those equations. 

A classification of second order integrable 2-component evolution equations 
has been given in |35j . The lemmas 6.3, 6.4, 6.5, and 6.6, proved in there, are 



special cases (n = 2) of Proposition 15 18 19 and 17 respectively. In pS] the 



full analysis of higher order symmetry conditions was carried out completely. 
From this it follows that all second order integrable equations with quadratic 
parts are derived from equations 2.3, 2.4, and 2.5. Note that the authors of [55] 
excluded equations of type 2.1 and 2.2. 

A classification of third order 2-component evolution equations with weight- 
ing (2, 2) and symmetries of order 5, 7, or 9, is given in [9 . Of the five equations 
listed [21 Theorem 3.3], two are nonlinear injective and have diagonalizable linear 
part: 

' \ _ / ^3 + UUl + VVl 

Vt J \ -2w3 - uvi 

and 

\ _ / 4u3 -I- 3w3 -t- Auui + vui + 2uvi 
Vt I \ 3u3 -f W3 — 2vvi — Avui — 2uvi 



(35) 



(36) 

When put in Jordan form the ratio of the coefficients of the linear part of equa- 



tion (36 1 becomes a/b — —30 — 2, where <j) is the golden ratio. Our diophantine 
approach perfectly explains the 'unusual' symmetry pattern. We remark that 
the conjugate of cj) gives rise to another equation with a/b = —3(1 — 0) — 2 = 
l/(— 30 — 2), which can be obtained by interchanging u and v. A similar remark 
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holds for all equations derived from the 'symmetric' Propositions 18 and 19 For 
example, by interchanging u and v in equation 4.11 we get an equation with 
a/6= 1/(7 + 47) = 7-47. 

According to Theorem 3.3] we have the following. A non-decouplable fifth 
order two component equation in the KDV weighting, possessing a generalized 
symmetry of order 7 can he reduced by a linear change of variables to a symmetry 
of lower order equations, or to the Zhou- Jiang -Jiang equation. 



Ut 
Vt 



( U5 — 5(2uw3 + 5itiU2) + 15(2wu3 + 3W1W2) + 20uiw^ \ 

— 30(miw^ + 2uvv\) 
— 9i'5 + h(2uzv + lu2V\ + 'du\V2 + 6UW3) — 10(2uuiw 

V \-2v?vi^'iv^vi) ) 



However, since the ratio of coefficients of the linear part ajh = —1/9 does not 
appear in our list, this equation has an approximate symmetry at lower order. In 
fact, the equation has a genuine symmetry at third order. The Zhou- Jiang- Jiang 
equation is in the hierarchy of the Drinfeld-Sokolov type [7] equation 



Ut \ _ / -3vvi 

Vt J \ V3~ uiv - 2uvi 



(37) 



which is linearly equivalent to a third order equation that appears in the same 
paper P Equation (17)], cf. [40, Sections 3.2.1, 4.2.6]. The special value of the 
ratio a/& = — 1/2 in equation ( |35| also does not appear in our list and is due to 
higher grading constraints. At the end of [24 two fifth order systems are given 
with ratios (9 - 5V3)(9 -I- ^VS)'^ = 26 - 15\/3 and -1/9. The first system 
derives from equation 5.6 and the latter from equation 5.3, with r a primitive 
6-th root of unity. 

We conclude with a more philosophical remark and some ideas on future 
research. The concept of generalized symmetry really is about local symmetry. 
The (inverse) Gel'fand-Dikii transformation translates polynomials in the sym- 
bols x,y into local differential functions, that is expressions in u,v and their 
derivatives. A question arises: can we also translate rational functions in the 
symbols x,y7 The answer is yes. One could think of non-local variables Ui,Vi 
with > J G Z. Here a negative index indicates integration and would be 
such that Dr^iui^i) — u^ for all i. We can expand rational functions in Taylor- 
series, which are transformed into non-local differential sums. For example, 
consider the rational function F = l/(xi -I- X2). Its symmetric Taylor-series is 

-, 00 t. JU 

fc=o ^2 ^1 

which is transformed into the non-local object 

00 

k=0 
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and we have D^F — v? . In this non-local setting every equation has a symmetry 
at any order. For example, the equation, with a ^ 1, 

Ut \ _ f o,ui + {I — a)uv 
Vt I \ vi + {1 — a)uv 

and its approximate symmetries (but the ones in A^''^ ® -4°'^), are in the ap- 
proximate hierarchy of the zeroth order non-local equation 

Ut \ _ ( u + UV^i 

Vt J \ V - u^iv 

Still, there would be a quest for equations, or symmetries, that are in certain 
sense close to local. 

Since the pioneering work of Sanders and Wang [34], apart from extending 
the result to equations with more more components [SJ |33 [H] , the symbolic 
method has been further developed in order to classify non-commutative [30] . 
non-evolutionary [5TJ|531[2S], non-local [1TJ[55], and multi-dimensional equations 
|1T] . In classifying non-local equations the concept of quasi-locahty, introduced 
in [25] is the key idea. Only certain types of non-localities are allowed by 
considering different extensions of the ring of differential polynomials. 'Nnon- 
evolutionary' equations are treated as evolution equations with possible non- 
localities. For example, a Boussinesq type scalar equation 

Utt = K{u, Ux,Uxx---,Ut,Uxt,--- 

can be represented as a two-component evolution equation, 

Ut \ ^ f V 
vt ) \ K{u,Ux, . . . ,v,Vx, . . .) 

However, different evolutionary representation may exist, and some might have 
local symmetries, whereas others might posses non-local symmetries |23j . An- 
other example is the Camassa-Holm type equation 

n^t \ ( cmUx + urrix 
Uxx J V u~m 

which is integrable when c = 2 ^ and c — 3 [6 . By eliminating the variable u 
this equation is written as 

TTif = cmAnix + m^Ani 

where A (1 - Dl)-^ is a non-local operator. Its symmetries are quasi-local 
expressions in D^-, and A-derivatives of m [2T1 124j . 

Another interesting problem would be to classify non-evolutionary equations 
as they are, that is, to apply the symbolic method for polynomials in both x and 
t derivatives, developed in |41j . In the setting of non-evolutionary equations, 
there is a clear distinction between the equation and its symmetries. When 
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A = represents the equation, a function Q is an infinitesimal symmetry if 
the prolongation of Q acting on A vanishes modulo A [23 Theorem 2.31]. 
For example, we have Utxx — UtUx as a symmetry of the Boussinesq equation 
Utt = ^xxxx ~ '^UxUxxj but not vise versa. Recently we formulated an equivalent 
symmetry condition in terms of Lie-brackets, in the setting of passive orthonomic 
systems [U]. This could play a role in the general classification problem. We 
also found that certain hierarchies of evolution equations appear as symmetries 
of non-evolutionary equations. For example, it seems that all symmetries of the 
Sawada-Kotera equation [36] 

5 3 

Ut = U5- SwgWi -I- -Wj^ 

are symmetries of Utxx = UtUx, which is a special case of Ito's equation |31j . 
Based on Lax-pairs a similar observation was made in jllj . 
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